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ABSTRACT 


The effect of delay in the control system input on the stability of 
a continuously acting controller which is designed without considering 
the delay is studied here. The stability analysis of a second order 
plant is studied analytically and verified numerically. For this example 
it is found that the system becomes unstable for a delay which is equiva- 
lent to only 16 percent of its natural period of motion. It is also ob- 
served that even a small amount of natural damping in the system can in- 
crease the amount of delay that can be tolerated before the onset of insta- 
bility. The delay problem is formulated in the discrete time domain and 
an analysis procedure suggested. The maximum principle from optimal control 
theory is applied to minimize the time required for the slewing of a general 
rigid spacecraft. The slewing motion need not be restricted to a single 
axis maneuver. The minimum slewing time is calculated based on a quasi- 
linearization algorithm for the resulting two point boundary value problem. 
Numerical examples based on the rigidized in-orbit model of the SCOLE also 
include the more general reflector line-of-sight slewing maneuvers. 
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I. INTRODUCTION 


The present grant extends the research effort initiated in previous 
grant years (May 1977 - Feb. 1986) and reported in Refs. 1-12 . Techni- 
ques for controlling both attitude and shape of very large inherently 
flexible proposed future spacecraft systems are being studied. Among 
possible proposed future applications of these large spacecraft systems 
(LSS) are: Earth observation and resource sensing systems; large 

scale multi-beam antenna systems (e.g. for use in mobile communications); 
orbitally based electronic mail transmission; and as in-orbit test models 
designed to compare the performance of flexible LSS systems with that pre- 
dicted based on computer simulations and/or scale model Earth-based labora- 
tory experiments. 

The present report is divided into two parts. This volume, designated 
as Part I,, contains four chapters. Chapter II is based on a paper pre- 
sented at the Sixth VPI§SU/ALAA Symposium on the Dynamics and Control of 
Large Structures which focuses on possible stability problems in LSS con- 
trol systems in the presence of delayed input which has not been taken into 
account in the design of the closed-loop system. For the special case of 
a hypothetical single degree of freedom system the effect of time delay is 
considered both analytically and numerically. The effect of inherent (natural) 
damping in the system is also analyzed. The control- problem with delayed 
input is also formulated in the discrete time domain. 

*References cited in this report are listed separately at the end of each 
chapter. 
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In the following chapter the problem of minimum time attitude slewing 
of a general rigid spacecraft is developed based on Pontryagin’s Maximum 
Principle. The slewing motion need not be restricted to slewing about 
a single axis, and the final attitude error can be made as small as required. 
The control torques and forces are computed and the minimum slewing time 
is determined using the quasilinearization algorithm for the resulting 
two -point -boundary value problem. Numerical examples based on the rigidized 
model of the Spacecraft Control Laboratory Experiment - SCOLE^ - 14 include 
both single axis slewing as well as the more general reflector -line -of 
sight slewing maneuvers . (A paper based on this chapter has been accepted 
for presentation at the ALAA 26^ .Aerospace Sciences Meeting, Reno, Nevada, 
Jan. 1988.) 

Finally Chapter TV describes the main general conclusions together 
with the future recommendations. The effort described here is being con- 
tinued during the 1987-88 grant year in accordance with our most recent 
proposal. 

It is planned that Part II of this report will be based on the Ph.D. 
dissertation entitled, "Ch the Dynamics and Control of the Spacecraft 
Control Laboratory Experiment (SCOLE) Class of Offset Flexible Systems," 
currently being prepared by Mr. Cheick Modibo Diarra. This document will 
focus on the modelling, stability analysis, and control law development of 
the SCOLE orbiting configuration, including the effects of the mast flexi- 
bility. 
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II. STABILITY ANALYSIS OF LARGE SPACE 
STRUCTURE CONTROL SYSTEMS WITH DELAYED INPUT 


Abstract 


large space structural systems, due to their inherent flexibility 
and low mass to area ratio, are represented by large dimensional mathe- 
matical models. For implementation of the control laws far such systems 
a finite a mo unt of time is required to evaluate the control signals; and 
this time delay may cause instability in the closed loop control system 
that was previously designed without taking the input delay into con- 
sideration. The stability analysis of a simple harmonic oscillator 
representing the equation of a single mode as a function of delay time 
is analyzed analytically and verified numerically. The effect of in- 
herent damping cn the delay is also analyzed. The control problem with 
delayed input is also formulated in the discrete time domain. 

I. Introduction 


Large flexible space structures have been proposed for possible use 
in communications , electronic orbital based mail systems, and solar 
energy collection. 2 The size and the low mass to area ratio of such 
systems warrant the consideration of the flexibility as the main contri- 
bution to the dynamics and control problem as compared to the inherently 
rigid nature of earlier spacecraft systems. For such large flexible sys- 
tems, both orientation and surface shape control may often be required. 

The equations of motion describing the shape of any large space 
structure are either represented by a few partial differential. equations 
or a large number of ordinary differential equations. As the partial 
differential equations are difficult to solve for control system design 
purposes, the structured dynamics are ccnmcnly described using Finite 
Element Methods (REM) . Two typical large space structures namely the 
Hoop/Coluirn antenna^ and the Space Station initial operational configura- 
tion (IOC) 4 are both described using 672 degrees of freedom. Thus the 
dynamics of a large: space structure can be written as^ : 

M z + K Z = U c (1) 
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where 


M = NXN mass/inertia syimetric matrix 

t 

K = NXN stiffness syimetric matrix 

Z — NXL generalize coordinates representing the degrees 
of freedom 


Uc* 


influence of the external forces in each degree of 
freedcm = B'U. 


T 

4> M4> 


With the modal transformation 

Z - q 

and the proper t ies of the modal transformation such as 

r 

<f> r K<j> * diag oj| 

and neglecting the higher modes, equation (1) can be written in standard 
state space form as 

X » AX + BD (2) 


where 


2nxl state vector representing modal coordinates 
and their velocities [q,q] T 


U = mxl 


oontrol vector 


A = 


0 

i -i 

ii 

1 nxn 

T 

i 


i 0 
2* 


—0) j 


J 


system matrix 


B = 


nxm 


<|> T B' 


nxm 


control influence matrix 


H. Control with Delayed Input 

The proposed control systems for large space structures are bgp e d 
on state variable feedback of the form: 


U = -FX 


( 3 ) 
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and the control gain matrix, F, is designed using techniques such as the 
linear quadratic regulator (LQR) theory®, pole placement', and/or linear 
quadratic Gaussian/loop transfer recovery (LQG/LTR) . 

For the case when the ..complete state is not available for feedback, 
an estimate of the state, X, is obtained using an appropriate estimator 
from the measurements of the form 

Y » CX (4) 


where 

Y - Hxl measurement vector 
C - JLxn sensor influence matrix 


In general, it is assumed that the estimated' state, X, is instanta- 
neously available.. As the state estimator is implemented using a 
digital computer and the number of the status (2n) is of the order of 
hundreds for a large space structure, the computational time becomes 
appreciable. Thus, in the present paper, the stability of the closed 
loop control system, with the control as given in equation (2), is 
analyzed as a function of the delay time (h) using the modified control 
law of the form: 


U(t) - -FX(t-h) (5) 

The characteristic equation of the closed loop system 

X - AX(t) -BFX(t-h) (6) 

is given by 


_qh 

G(s,h) - det (sI-A+BFe ) - 0 
which, in turn, can be written as 

G(s,h) » 2 £ P i (s)e" shi * 0. 

i-0 


(7) 


( 8 ) 


The roots of the characteristic equation, (8), as a function of the delay, 
h, are obtained from the corresponding auxiliary equation^ 


where 


G’ (s,h) - Z P. (s)(l-Ts) 21 (1+Ts) 4n “ 2i - 0 
i-0 

-sh _ r 1-sT -1 2 
e *• 1+sT ■* . 


(9) 

( 10 ) 
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Hie value of T for which the roots of the equation (9) cross the 
imaginary axis in the s-plane is obtained and the correspo nd! h is 
evaluated using the relation, (10). 


III. Example of a Fanncnic Oscillator 

The equation of motion representing the ith structural mode is 
the familiar harmonic oscillator and is given by 


x i + u i x i 


Considering the d el ayed velocity feedback of the form 


f ^ (t-h) 


with 


* 6/ * 0.5, 

the characteristic equation is given by 
G(s,h) = s 2 +36+6se“ sh = 0 


"iio V 3 ’® -3 * 11 = 0 


where 


P„(s) = s 2 + 36 


P x (s) = 6s 

The corresponding auxilary eq uati on is given by 


( 11 ) 


( 12 ) 


(13) 


£ P, (s) (l-Ts) 21 (1-rtS) 2 " 21 =0 (14) 

i=0 1 

1 ‘ e * (s 2 +36) (1-HTs) 2 +6s (l^Ts) 2 = 0 

or T^s 4 + (2T + 6T 2 ) s 3 +(l+36T 2 -12T)s 2 

+ (72Tt6)s+36 =0 (15) 

Using the Routh-Hurwitz criterion, it can be found that the roots of 
equation (15) cross the imaginary axis at u> = 9.7 for T = 0.0426. 

The corresponding dela y (h) can be calculated from the relation (10) 
with s = jto and is 0.16. This result can also be verified directly 
for this sinple system with the substitution s=jaj into eq uatio n (13) 1^, 
res u lt i ng in the value of oj and delay h for which the roots of the 
characteristic equation cross the imaginary axis. 
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Thus, equation (13) can be written as (keeping ^ and oj^) : 
(u^-ta 2 ) + j (2? jO)^ai)e ^ Co * 1 = 0 


2 2 , 


(u^ -a) + 25 ^u^a) sinuh) + j2£^u^u cosuh ■ 0 
For equation (17) to be satisfied 
cos uh ■ 0 or uh » 11 / £ 


and 


or 


2 2 

u^ -a) +2c^UjU ■ 0 

u ■ £± u ± ± “j/ 1+5 


Taking the positive value for u, the delay h, is given by 

ir/ 0 


<^[ 5 *+ / 1 + 5 ^ ] 


(16) 

(17) 

(18) 
(19) 


( 20 ) 


The value of h for 5 ± - 0.5 and u ± - 6 is 0.16 and thus the earlier 
result is verified. It is observed that an increase in damping reduces 
the tolerable delay (h) in the input. 

The equation of motion of a single mode with inherent (natural) 
damping and velocity feedback can be written as: 


X+2?^u) i X-h«j 1 X « f 


-2c 1 o» i X(t-h) 


( 21 ) 


where 5^ is the inherent damping ratio. 

The corresponding characteristic equation is given by 


• s 2 +25[d) i s-hD i 2 +2c i u) i 3e" sh 


0 . 


After substituting s * ju, equation (22) can be written as: 

2 2 t 

(<ik -ui +2 c^oku> sinuh) + j (25^u^u+25^u^u'cosuh) =* C 
For equation (23) to be satisfied for all u and h, we have 

2 ^i t0 i +2j *i tll i cosa) k 3 0 
or 

cos uh - -Cj/ 


( 22 ) 

(23) 

(24) 

(25) 
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Thus, for cosuih * < 1, the inherent damping must be less than damping due 
to control for instability. For £^<2;^, the system will always be stable. 

With the value of wh from equation (25) the frequency ui can be 
calculated, as: 

<* * m . O 5 2 - S '2 + ✓ 1+^2 - 5 '2 ] ( 26 ) 

i i i i 

and selecting the positive value of w, h is given by: 
cos 1 (“?i/C 1 ) ‘ 

h - 


01,0 sk' 2 + / i+c‘ic ,a 3 (27) 

ii i i 

For Zj * it can be seen that the delay, h, is half the undamped 
natural period of vibration. As the damping due to control increases, 
the tolerable delay (h) decreases and is in accordance with the observa- 
tion made in the case without the inherent damping. The effect of inher- 
ent damping in the system is to increase the amount of delay that the 
system can tolerate without become unstable as compared to the case 
without damping 

17. Discrete Time Domain 

As the controller is implemented on a digital computer, it may be 
more natural to consider the delayed input problem in the discrete time 
domain. 

The equations of motion as given by equation (2) can be written in 
the discrete time domain as 


where 


X(i+1) - A d X(i)+B d U(i) 


M 

A 


e M , B d - f 
U 0 

discretization time. 


B dt 


(28) 


The delayed input. problem can be considered in discrete time in one of 
the two following ways: 

i) Designing the controller of the form U = -FX('i) without taking 
into consideration the delay and then examining the effect 
of delay on the stability of the closed-loop control system. 

The control gain matrix F is designed such that the matrix (A^-B^F) 
has the eigenvalues within the unit circle. Then the delay is introduced 
into the control law as: 
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and 


U(i) - - FX(i-£) 

X(i+1) - A d X(i)-B d FX(i-H) . 


(29) 


(30) 


The stability of equation (30) can be studied; using the augmented system 
given by 




' - 

. 1 

1 

X(i+1) 


A d 0 0 0 -B d F 


X(i) 

X(i) 


I 0 0 0 0 


X(i-l) 

x(i-a+D 


0 0 0 I 0 


X(i-i)_ 



u 

A . 

Z(i-Kl) 


1 

w . 

Z(i) 


(31) 


or 


Z(i+1) - A d Z(i) (25) 

(ii) Designing the control by taking into account the delay in 
the input. 6 * 11 

Equation (28) can be modified as : 

X(i+1) -■ A d X(i) + B d U(i-£) (32) 

f\, 

The control law of the form U(i) * -FZ(i) can be designed from the 
augmented system: 


U(i) 


X(i+1) 


A, 0 0 o B J 

a a 


’ X(i) 


0 

U(i) 

- 

0 0 0 0 0 


U(i-l) 

+ 

I 

U(i-l) 


0 10 0 0 




0 

U(i-i+l)_ 


r 

o 

w 

o 

o 

o 

1 


_ U(i-E). 


0 


(33) 


<v 

Z(i) 


or 


Z(i+1) - A d Z(i) + B d U(i) 
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Thus the input U(i-i) is a function of the previous inputs, U(i-i-l), 
U(i-&-2),..., and the previous states X(i-Jl). Though this design can 
take delay into consideration, the sequence of the control signals: 
U(i-i), U(i-JM-l) , must be generated at an interval of one step and, 
thus, the original delay problem is not completely solved. 


Conclusions 


The effect of delay in the input on the stability of the continuous 
time controller that is designed without taking this delay into consider- 
ation is presented. The closed-loop control system of a second order 
plant becomes unstable for a delay of 0.16 seconds, which is only 16 per- 
cent of its natural period of motion. It is also observed that even a 
small amount of inherent (natural) damping in the system can increase the 
amount of delay that can be tolerated without the system becoming unsta- 
ble. The delay problem is formulated in the discrete time domain and an 
analysis procedure is suggested. 
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Appendix - Chapter II 

Stability Analysis of Second Order System 
with Delayed State Feedback 

% 

As a second order differential equation describes the dynamics of 
a single mode of any large space structure, the stability analysis of such 
a system with delayed state feedback is analyzed and the amount of delay 
that can be tolerated by the system without becoming unstable is arrived 
at analytically. 

The differential equation of second order with state feedback can be 
written as: 

x i + 2?ju> i x i + aj i 2 x i = -k r x i (t-h) - kpX^t-h) (1) 

where 

x^ = i^ 1 modal coordinate 
Wj, = i^ 1 natural frequency 

= itk mode inherent damping ratio 
k r = rate feedback gain 
kp = position feedback gain 
h = time delay 

The feedback gains k r , kp are designed for the required stability and 
transient response specifications without taking the delay into consideration. 

The inherent damping ratio, q' and the feedback gains, k r and kpj will 
give rise to five possible combinations as shown in Table 1 and are thus 
analyzed separately for mathematical convenience and easy understanding. 
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Case I: 


q = 0, kp = 0 and- k r > 0 
The differential equation of the system can be written as: 

x i + “i x i = “ k r ^if t_h ) ( 3 ) 


Case 

Si 

k r 

k 

P 

I 

= 0 

> 0 

= 0 

II 

> 0 

> 0 

= 0 

III 

= 0 

> 0 

> 0 

IV 

> 0 

= 0 

t 0 

V 

> 0 

> 0 

+ 0 


Note: The remaining three combinations are 

neither feasible nor of interest. 


Table 1: Feasible Combinations of q , k r , k_ 

for Stability Analysis P 

and the corresponding characters tic equation is given by: 

s 2 + a 2 + Zqa^se" 5 * 1 = 0 C4a) 

where k r = 2s^. 

The value of h for which the roots of equation (3) cross the imaginary 
axis can be evaluated by substituting s = ju). 

Thus 2 

uk -u> + j2qui^u)Sina>h + 2s^uKu)coswh - 0 (4b) 

For equation (4b) to be satisfied 

sinwh = 0 
2 2 

and u>. -uj + 2t-u)-w cosuh = 0 


( 5 ) 



Thus 


and 


toh = V 2 
h = 


*/2 


^ [? i + ^ ] 

Case II: > 0, k r = 2?^^ and kp = 0 

The characteristic equation of the system described by equation (1) 


( 6 ) 


is given by 
2 2 


(co^-oj +2?^aj^(jjsinu)h) + j (2s|u)^w+2c^aj^(ocoS(i)h) = 0 


(7) 


Thus 

and 


coswh = 


( 8 ) 


^ _ cos~l(£j/£i 

»i [ /If^rr + / 3 

For the case where thesystem will always be stable since no value 

of h exists for which the roots of (7) cross the imaginary axis. A plot 
of ujh versus for various values of ^ is shown in Figure 2.1. 

Case III: c! » 0, kp = k r > 0 

The characteristic equation is given by 

s2+u)2+k r se' s * 1 +kpe" s k = o (9) 

(u)?-u)?-hok r sinuh+kpCOSa)h) + j ( a»k r cosu»h-kpSinajh) = 0 (10) 

tan <uh = 


or 

Thus 


tok-j* 


and 


O) 2 = ^[(2u>|+k2) + /Ck^+4aj2k^+4k^] ] (11) 

Plots of hw^ versus k r /cu.^ for various values of kp/oi? are shown in 
Figure 2.2, It can be seen here that these are many combinations of kp and 
k r for which the roots of Eq. (10) can cross the imaginary axis - i.e. 
value of hu^ which leads to instability. 
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Case IV: 


> , k r = 0, kp ? 0 
The characteristic equation is given by 
(u>^ -co^+kpCoscoh) + j (Zt^'aj^ui-kpSinwh) = 0 


Thus 

and 


sinwh = 


2q co^oj 


■¥ 


2 = 2 ri _ ? r , 2 ^ A 2 y rri _ 0 , f 2 ~ 2 - 


= a, i (l-2? i ) + «>/[(l-2ci Z ) + CkpA-i 2 ) 1 
The plots of hu^ versus kp/co^ f or various values of are shown in 
Figure 2.3 


( 12 ) 

(13) 

(14) 


Case V: q > 0, k t > 0, kp f p 

The characteristic equation is given by 

7 7 

(ui^-to +o>k r sinuh+kpCOSwh) 

+ j(2q a) +oik r cosojh-kps inoih) = 0 
By equating the imaginary part to zero, uh can be evaluated as 


> u • -1 2?- 

um = sin x c i 

t "VI "V 


U> .ID 


T 


kp +oj k r z 


-1 r Mk r . 
- tan • L ( — t — ) 


(15) 


(16) 


y 

after substituting wh in the real part of equation (15) > 

the following equation in the single unknown variable co can be obtained 


i^-u^+wk r sin (sin”V - tan ~^( Mkr )) 

kp 

-1 _1 Ojk-r 

kp cos (sin y- tsn ( — ^ — )) = 0 


(17) 
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Using equations (17) and (16), the limiting value for given values of ' , 
kj., kp and u>^ can be determined. As the equation (17) is nonlinear, 
numerical procedures may have to be used and thus the generalized plots 
similar to the other cases may be obtained; 
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0,000 * 150 ,300 *450 .600 .750 .900 1.050 1 .200 1.350 

active damping ratio 5 

Figure 2.1: Plots of hoi- vs correspondence to Case II with . 5 ! 

as a parameter. L 
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III. minimum time attitude slewing 

MANEUVER OF A RIGID SPACECRAFT 

1 . Introduction 

The problems of large-angle attitude maneuvers of a spacecraft 
have gained much consideration in recent years E 1—8 , 11], i n these 
papers, the configurations of the spacecraft considered are: (1) com- 
pletely rigid, (2) a combination of rigid and flexible parts, or (3) 
gyrostat-type systems. The performance indices usually include minimum 
torque integration, power criterion, and frequency-shaped cost func- 
tionals, etc. Also some of these papers used feedback control tech- 
niques. In this paper, we try to concentrate on the minimum time slew- 
ing problem of a rigid spacecraft. 

In Ref. [2], the author studied the rapid torque-limited slewing 
of SCOLE about a single axis (x-axis) about which the spacecraft has 
a small moment of inertia. The control torque about this axis is of a 

bang-bang type or a bang-pause-bang type. The author computed the 

* 

slewing motion on the simplified model of the rigidized SCOLE EH, 
then worked on the practical rigidized model (with nonzero products 
of inertia); hence, this leads to a large error of the attitude 
after the slewing. Also it seems that no details were given for the 
controls about the other two axes (y, z). 

In the present paper, we apply optimal control theory (Maximum 
Principle) to the slewing motion of a general rigid spacecraft (in- 
clude the rigidized SCOLE, without simplification). The slewing 
motion need not be restricted to a single-axis slewing. The attitude 
error at the end of the slewing can be made as small as required. 

All the controls (torques and forces) are computed and the minimum 
slewing time is found by using the quasilinearization algorithm for 
the resulting two-point-boundary-value problem. 
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2. Attitude Description and State Equations 


2.1 Attitude Description and Euler Rotation 

Let a * [ai a 2 £ 13 ] T reprensent a set of unit, orthogonal vectors 
of an inertial reference system, and b = [b^ b 2 b 3 ] T a set of unit, 
orthogonal vectors of a body-fixed coordinate system of a spacecraft. 
Then, the attitude of the spacecraft relative to a can be described 
by a direction cosine matrix C with C satisfing the relation 


b = C a 


( 1 ) 


and 


C 


2 2 2 2 
^0 +q l- q 2“ q 3 

2 < q l q 2 _q 0 q 3 ) 

2(q 1 < 3 3 +qo < 3 2 ) 


2(q 1 q 2 +q o q 3) 

2 2 2 2 

q 0 +q 2“ q 3“ q l 

2(q 2 q 3 - < 3o q i) 


2(q 1 q 3 -q 0 q 2 ) 

2 ( q 2 q 3 +q 0 q l > 
2 2 2 2 
q 0 +q 3“ q l“ q 2 


(2) 


where q=[qo qi q2 q 3 l T is the attitude quaternion vector and subject 
to a constraint equation 

q T q=l ( 3 ) 


It can be seen that q can be used not only to represent an atti- 
tude orientation of a spacecraft, but also to describe a rotation of a 
rigid body (spacecraft). For example, when a rigid spacecraft rotates 
about an axis defined by a unit vector £ = [£3 £2 ^ 3 ^ T fixed in both 
a and b, the quaternion describing this rotation is 
' qr = cos ( 0 / 2 ) 

(4) 

qi = £i sin (®/2) i— 1 ,2,3 


where 0 is the rotation angle. 
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The Euler rotation theorem tells us that an arbitrary orientation 


of a rigid body can be accomplished by rotating it about a certain 
eigenaxis, £ = [£^ &2 through 6 angle from its initial position* 

By means of this theorem we can find the desired rotation quaternion, 
q, between the initial position q(0) and the final orientation q(tf) 
by the relation 


(5) 


■ ^0 ' 


q oo q 10 q 2 o S30 


1 

m 

0 

O’ 

1 

Si 


■ q io q 00 q 30 -q 2 o 


q lf 

S 2 


“ q 20 “ q 30 q oo q io 


q 2f 

.S3 


. “ q 30 q 20 - q 10 q 00 . 


. q 3f . 


where the second subscript "q" and "f" represent the initial time 
and final time, respectively. From here, we can find £ and 0 


0=2 arc cos qg 

£i “ qi/^l-V 


(6) 


i-1,2,3 
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The attitude quaternion and the angular velocity of a rigid 
spacecraft satisfy the following kinematical and Euler dynamical 


equations . 


q = ( 1 / 2 ) 00 q 
I = CO I CL> + B u 


* (7) 
( 8 ) 


where 


U) angular velocity vector in the body system, = (*) 2 0)3 ] T 

u control torque and force vector, u = [ u ^ u 2 u 3 ... u n ] T 


0 -(02 -(0 3 
( 0 i 0 U ) 3 -U^ 

U) 2 -(dg 0 
CO 3 ^2 ”^L 0 


'0 (O3 - 60 ? 

UO - -0)3 0 U>i , 

U?2 0 


I - -I 


11 

- I 12 

1 

H 

I-* 

OJ 

12 

I 22 

_I 23 

13 

* T 

" 23 

X 33 


and B is a 3xn alignment matrix. Eq.( 8 ) can be rewritten, by pre- 
multiplying the inverse of I, as 

00 = I"1 CO I uO + I"1 B u (9) 

The associated initial and terminal boundary conditions of the 
states, q, ( 0 , are prescribed: 

q ( t =0 ) , C 0 (t= 0 ); ( 10 a) 

q(t f ) , 0 )(tf) ( 10 b) 
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3. Optimal Control — Two Point Boundary Value Problem 
In this paper, we try to minimize the slewing time tf, 

f tf 

tf - \ dt (11) 

' 0 

under the constraints that the elements of the control vector u have 
their upper and lower limits, respectively 

u jmin u j £. u jmax * j = l,2,3,...,n (12) 

Generally speaking, minimization of tf under the constraints 
(12) will result in a so-called two point boundary value problem in 
which several controls (at least one ) will reach their bounds during 
the slewing time, tf. To explain this point, let us first consider a 
well-known special case where there are only 3 control torques, u i, 
u 2' u 3» about the 3 principal axes of the spacecraft, respectively 
(i.e. diagonal matrix I ). For this case the minimum time rotation of 
the spacecraft about one of its principal axes will yield the 
following results: the control torque about this axis is of a 
bang-bang type, while the other two torques remain zero. Otherwise, 
if the slewing motion is not about a principal axis, none of the 3 
controls remain zero, but we can reason that at least one of the 3 
control inputs reaches its bounds, except some jumps at the switching 
points during the period, tf. As for a general case where the control 
torques are about a body axis system which does not coincide with the 
principal axes (non-diagonal I) and some additional control forces, 
u 4' u 5» •••» u n» are available, the control laws become more 
complicated. 
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To handle the problem in which some controls reach their bounds 
and others do not, we introduce an additional cost function 

tf 


= j r f u* 

2 JO 


Ru dt 


(13) 


where u is the control vector, R is a proper weighting matrix. From 
Refs. [3] and [8], we can see that, for the case of rest-to-rest 
(i.e. (aXOJ^O, id(tf)=0) slewing with only 3 control inputs involved, 
if we use only (13) as a criterion and tf is long enough, the control 
torques are approximately linear functions of time, and the controls 
will not reach their saturation levels. But if we shorten tf in order 
to find a minimum time, some of the controls must reach their bounds 
and, thus, contribute more effort to the slewing. By continuing the 
shortening of tf, we can get a particular value, t*, during which 
at least one of the controls remains as bang-bang with one switching 
point, while others are generally not of the bang-bang type. This 
value, t*, is called the minimum time which is required. 

The motivation for using (13) as our cost function are: 

1) Ease of using the quasilinearization algorithm. 

2) No need to determine the switching points. 

3) Easy to guess the unknown initial values of the costates. 

3.1 Necessary Conditions 

The Hamiltonian, H, for the system (7), (9) and (13) is 
H = (1/2)u t Ru + p T q + r T iO 

= (1/2)u t Ru + (l/2)p T 3q + r T (I _1 & I (0 + I' 1 B u) (14) 
where p and r are costate vectors associated with q and i0. 


P = [P 0 Pl P2 p 3 1 T - 


r - [ r ! r 2 r 3 ) T - 
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By means of the maximum principle, the necessary conditions for 
minimizing J, are 


p = - Qh/ 2 q) , ==> 

p - (1/2) % p 

(15) 

r = - Oh/^aJ) , ==> 

r = g(W, r) + (l/2)[q]p 

(16) 


where g(C0, r) is a 3x1 vector function of ^ and r, and the detailed 
form of g(W, r) can be found in Appendix I; [q] is a 3x4 matrix 



q l 


-q 3 

q 2 

[ q ] - 

q 2 

q 3 

■ q o 

-q x 


. q 3 

“ q 2 

q i 

- q o 


The initial values of p, r are unknown, p(t=0), r(t=0). 

If u is a 3x1 control torque vector and B is a 3x3 nonsingular 
matrix, R can be a positive-definite matrix defined by 

' R = B T B (17) 

R u +B t I -1 r = 0 
u = - R -1 B t I" 1 r 

= - B' 1 I“1 r (18) 

if u j < u j m i n ; 

if u jmin 1 u j i u jmax > (19) 

if u j > u jmax* 
j = l f 2 , 3 . 

If u is an nxl (n>3) vector, B is a 3xn matrix, the R formed by 
(17) is a semi -positive-definite matrix. To circumvent the singularity 
of R, we introduce a 3x1 vector, v, 

v = B u (20) 




3H 

du 


■ 0 


or 


we have 


f u . . 

r jmm • 

u j = < -(B" 1 I" 1 r ) j , 


u . 

jmax» 
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Then 


(2h/$v) = 0 , ==> v = -I -1 r (21) 

By means of pseudo-inverse of matrix B, B + , we can get u 

u = B + v = B t (BB t ) _1 v 

= - B^BB^-iCE-ir) (22) 

The control laws are 



f u . . 

[ ]mm ' 

if u. < u . . . 
j v }min' 


u j , . 

f - (B + I -1 r)j. 

if u. . rU.rU. 

11 jmm j i jmax • 

z (23) 


1 u • 

\ }max' 

if u . > u . 

j jmax* 




3 1 # 2 / • • • f n • 



Note that Eq.(23) is reduced to Eq.(19) if B“1 exists. 

In summary, we seek the function q(t), tt(t), u(t), p(t), and r(t) 
which satify the equations (7), (9), (15-16), (23) subject to the 
boundary contions (10). 

3.2 Properties of the initial values of p 

The key to settle this problem is to find the unknown initial 
values of the costates 

p(0)=[ Pqo p 10 P 20 P 30 ] T and r(0)=[ r xo r 20 r 3()l T 
Notice that the coefficient matrix of the right side of Eq.(15) is 
anti-symmetric, so, 

p T p = 0 

i.e. p T p = constant 

The extra constant is usually treated as an unknown and is determined 
by iteration. This results in more computational effort. However, as 
we shall prove, this unknown constant can be easily selected without 
changes in the optimal controls [83. 
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Compare Eqs.(7) and (15), they have the same coefficient matrix 
on the right sides. Therefore, they have the same state transition 
matrix. Let Q represent this 4x4 matrix, then the q and p at any 
instant can be obtained by 

q = Q q( 0 ) , p - Q p(0) (24) 

We know that Q satisfies the following matrix differential equation 


Q = (1/2) 0) Q 

Ref. [10] shows that Q, the solution of (25), has the form 


(25) 






q ll 

- q 12 

- q 13 

“ q 14 " 










q 12 

q ll 

q 14 

- q 13 








Q - 


q 13 

' q 14 

q ll 

q 12 





(26) 





L q 14 

q 13 

- q 12 

q ll. 






On 

substituting Eq 

;.(26) into Eq. 

(25) we 

s can verify that only 4 of 

the 

16 

are independent. We 

rewrite the 

first 

equation 

of (24) as 




f* % 


■* 











q 0 


q oo 

- q io 

“ q 20 

_q 30 


q ll 






q l 


q io 

q oo 

- q 30 

q 20 


q 12 



(27) 



q 2 


q 20 

q 30 

q oo 

- q io 


q 13 

• 





. q 3. 


. q 30 

_q 20 

q io 

q oo ' 


. q 14. 




where 4oo aq o(°)> q 

[ 10“ q l(° ) ' 

etc. 

It is 

clear 

that the 

coefficient 


matrix of Eq.(27) 

is orthogonal, 

so 











r _ 




1 

r _ 





' q il 


q oo 

q io 

q 20 

q 30 

q 12 


" q 10 

q oo 

q 30 

- q 20 

q 13 


_q 20 

- q 30 

q oo 

q io 

. q 14. 


. " q 30 

q 20 

“ q 10 

q oo 


*0 

*1 


l q 3 


(28) 
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From Eq.(28) we get 


q ll + q 12 + q 13 + q 14 = ! 

This means that Q is also orthogonal. On the other hand, we have a 

similar equation for p, 

P 0 P 00 " p 10 “ p 20 “ p 30 q ll 

P 1 p 10 p 00 ” p 30 p 20 q l2 

■ (29) 

P 2 P 20 P 30 P 00 “ P 10 q 13 

P 3, . P 30 “ P 20 p 10 P 00. . q 14. 

After substituting Eq.(28) into Eq.(29) and eliminating <3n, <3 12 , 

<313, ^14, one arrives at 




associated costates. Eq.(30) can be rewritten as 



( 32 ) 
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Substituting of Eq.(32) into Eq. (16) results in 
r = g(U, r) + 



or 


r * g(td, r) - 



0 

2 2 2 2 



V 


0 

^1^2- q 3 

2( q l q 2 +q 0 q 3^ 

2 ( q l q 3~ q 0 q 2 ) 



1 



2 2 2 2 



d l 


0 

2 ( q l q 2“ q 0 q 3 ) 

q 0 +q 2” q 3'" q l 

2(q 2 g 3 +q o q i) 



2 




2 2 2 2 


d2 


0 

2( q i q 3+ q o q 2) 

2 ( q 2 q 3~ q 0 q l ) 

q 0 +q 3 -q l“ q 2 








d3 


or 

r - g (l*) t r) - (1/2) C d (33) 

where d=[d^ d 2 ^3] T / c is just the attitude matrix given by Eq.(2). 

It can be seen that r is independent of from Eq. (33), and u 

depends only on r, from Eq. (23). Therefore, u is also independent of 
^ 0 - This means the arbitrary selection of the value of yields the 
same extremum control, u. Now- we can explain the results in Ref. [11]. 
In view of Eq.(31), we have 


jZ ,z .z z z z z 

+ d i + d 2 + d 3 = p Q0 + p 10 + p 20 + P 3Q 


(34) 


If we set ^0“® the norm of the initial costates in Eq.(34) reaches 
its minimum, the solution of which is considered in Ref. [3]. From 
Eq. (31) we can also know that ^o = 0 means 

p(0)Tq(0)=0 
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4. Initial Values of Costates and the Slewing Time 

By means of Euler's eigenaxis rotation theorem, from the known 
attitudes at the initial and final time, q ( 0 ) and q(tf), we can find 
a unit vector (eigenaxis), £ , which is fixed in both the body axes and 
inertial coordinate system, and a rotation angle, 0*. Then the atti- 
tude changes from q(0) to q(tf) can be realized by rotating the 
spacecraft about the axis, £, through the angle, 0*. 

Theoretically, there are many ways through which we can change 
the attitude from its initial value, q(0), to its final value, q(tf). 
For example, this change of attitude can be achieved by successively 
rotating the spacecraft about the x, y, z axes (i.e. 1-2-3 rotations) 
through certain displacements in the angles, ®3» respectively. 

To do this way, we need to speed up (and slow down) the spacecraft 3 
times, and the total rotation angle is, e i+ ®2 + ®3* On the other hand, 
for the Euler rotation, we only need to rotate the spacecraft about £ 
once through the angle 0* which is less than the total angle required 
by any other way. Since the Euler rotation is simple and requires a 
smaller angle, it may take less time and consume less energy (torques 
and forces). Therefore, in view of our cost functions, (11) and (13), 
it is reasonable to think that the optimal slewing is near the Euler 
rotation. We shall call this rotation the "expected rotation", which 
is determined only from the initial and final attitude of the space- 
craft and will be used in obtaining a set of approximate unknown 
initial values of the costates and the starting solution of the 
quasilinearization algorithm. 

4.1 Initial Values of Costates 

Before starting the quasilinearization algorithm, we need to 


3-12 



guess the unknown initial values of the costates, p and r. Considering 
the analytical solution about a single principal axis maneuver in Ref. 
[3], we define a rotation angle 0(t), about an arbitrary axis t, 


e(t)=e(o)+ e(0)t + - e(o)t 2 + - e(0)t 3 

2 6 


(35) 


where 0(0), -0(0), 0(0), 0(0) are constants to be determined. 

For simplicity, here we only consider the solution of © with the 
following boundary conditions 

0(O)=O, 0(O)*O, 0(t f )=0* 0(t f )=0 (36) 

These conditions correspond to the boundary conditions of the states 
q( 0 ) , U(0)*0, q(t f ) , <d(t f )=0 

Substituting Eq.(36) into Eq.(35) yields 

0(0)= (60*/tf) - (40(O)/t f ) (37a) 

A 

3 2 

0(0)= - ( 120*/t f ) + ( 60(0 )/t f ) (37b) 

For the Euler's rotation, the angular velocity and its deriva- 
tives are expressed as follows 

CO = g©, LO = g©‘, LO = £’© (38) 

To approximately determine the initial values of p and r, we 
need to use the dynamical Eqs.(9) and (33). Upon using Eqs.(20) and 
(21), substituting v (u) into Eq. (9) and solving for r, we get 

r = I <0 ICJ - I 2 (39) 

and the derivatives 

^ = ^— ( i cO i co) - 1 2 60 (40) 

dt 

At the same time, from Eq.(33), and noting that C T =C“1, 

d = 2 C T [ g(0), r) - r ] (41) 
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At the time t-0, by putting Eqs. (37-38) into Eqs. (39-41), we can get 
the approximate values of r(0) and d. As for p(0), we can set Po(0)=0 
(since can be arbitrarily chosen) and solve for and Pf(0), i=l , 
2,3, by using Eq.(32). 

Now we determine 9(0) from known initial value Cd(0). Generally, 
CU(0) is not equal to £0 because £ is independent of 6J(0). Let e be 
the difference between them 

e = £0(0) - CJ(0) 

To find a minimum value of e T e, we differentiate e T e with respect to 
0(0) and note that £ T £=1 , we get 

0(0) = £ T (0(0) (42) 

By using the initial values p(0) and r(0) obtained above, and 
integrating the differential equations (7), (9), (15-16), with Bu in 
Eq . ( 9 ) replaced by v in Eqs. (20) and (21), we can get a set of values, 
q(t), <d(t), p(t), and r(t), 0 < t < tf, which will be used as the 
starting values of the quasilinearization method. 

4.2 Initial Value of tf 

The starting value tf(°) needs to be # made as close to the minimum 
time, tf*, as possible. This can be done by using the techniques simi- 
lar to those described above. Suppose the slewing motion is an Euler 
rotation about a vector, £, through an angle, 0(t). Then, by putting 
the first two equations of Eq.(38) and Eq.(20) into Eq.(8), we get 

l£0 = 0 2 £ , l£+v (43) 

For simplicity, we only consider the case v imin =_v imax* Then, let 
c i =v imax an d v i= c fT* i * the above vector equation can be expressed as 
the following 3 similar equations for 0(t): 


3 - 14 - 


i =1,2,3 


(44) 


0 fc>i Cjl C i 

where a^, b^, and are constants, -c^ is the normalized control 
about the ith body axis and 

| Ti | <1 i— 1,2,3 (45) 

Each equation of Eq.(44) with the boundary condition Eq. (36). can be 
treated as a minimum time control problem with the constraint (45). 

It is easy to see that the control for this problem is of a bang-bang 
type and the problem can be solved analytically to get the minimum 
time t*fi (i=l, 2, 3) as functions of 0* and 0(0). The results are 
shown in Appendix II. 

Since the only minimum time, tf*, that every equation of Eq.(44) 
can accept at the same time is the longest one, we use the largest 
one as our initial guess for tf. 

4.3 The Quasilinearization Algorithm ^ 9 ] 

We choose the quasilinearization algorithm to solve the two 
point boundary value problem because this method needs only to solve 
linear differential equations and it converges quadratically. 

Let z=[qo qi q2 93 W 1 w 2 w 3 P0 Pi P2 P3 r l r 2 r 3l T ' then Eqs.(7), 
(9), (15-16) can be replaced by 


z = f (z,u) 


(46) 

with the boundary conditions 



Zf (0 ) = Lf 

i=l,2, ***,8 

(47a) 

Zj(tf) * Mj 

j=l » 2 , •••,7 

(47b) 


The 8th initial condition is 2 q(0)=Po( 0) (arbitrarily chosen). 


Suppose that at the Nth stage of the iteration, an approximate 
solution reasonably close to the exact solution has been obtained. 
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Then, the linearized form of Eq.(46) at the (N+l)st stage of the 
iteration about the nominal functions, which are the solution obtained 
at the Nth stage of the iteration, is 

z = Of/dz)(z<N), u (N)) ( z (N+l ) _ z (N) ) 

+ (df/^u) ( z (N) ,u^ N ^ ) (u( N+ D-u( N ) ) + f(z( N ),u( N M (48) 
where ^f/dz and ^f/^u are the gradient matrices evaluated at z( N ), and 
u<»>. The u<N+l) in Eq.(48) is determined as follows. At a point on 
the optimal trajectory, it is assumed that u^ N ) satifies Eq.(23). If 
for some j (j-1,2, n), u^^j^uj^ (or Uj m i n ), then 

U (N+1). = U (N). (49) 

On the other hand, the linearized form of Eq.(23) is 

u (N+l)_ u (N)sa - B + I-l(r( N+1 )-r( N ) ) (50) 

i For i¥j (i=l,2, • * * ,n ) , 

I 

' u( N )jL - -(B + I~ 1 r< N ))i 

By putting this equation into Eq.(50), we have 
j u (N+l) i » _( B + x -l r (N+l) ) ± (51) 

Replacing u^ N+ D in Eq.(48) by Eq.(49) and Eq.(51), we can get a 
linear differential equation for z( N+ D 

z (N+l ) „ a(2 (N), u (N) )z (N+1) + d(z(»),u<H)) (52) 

which is nonhomogeneous and subject to the boundary conditions (47): 

z< H >i(0) = Zi (0) = L if i-1,2, *••, 8 (53a) 

z (N)j(t f ) = Zj (tf) =Mj, j-1,2, •••, 7 (53b) 

and can be solved using the method of particular solutions. 

Let yi, Y2' • ••» Y6 represent the solutions of Eq.(52) with the 
following initial conditions 


3-16 



yi(0)=[Li L 2 *** Lg 1 0 0 ••* 0] T 

Y 2 ( 0 ) = [Li 1*2 *'• Lg 0 1 0 •*• 0] T 

• • 

• • 

Y6(0)=[Li L 2 **• Lg 0 0 ••• 0 1] T 

Let y 7 be the solution of Eq.(52) with the initial condition 

Y 7 ( 0 )=[Li L 2 *•* Lg 0 0 ••• 0 0 ]T 

We now have 7 particular solutions of Eq.(52), Yi(t), y 2 (t), * • • , 
Y 7 (t), 0 < t < tf. The general solution of Eq.(52) can be written as 
a linear combination of these particular solutions 

z = BxYl +B 2 Y 2 + **• +B 6 Y 6 +B 7Y7 (54) 

where B=[Bx B 2 •** B 7 ] T are constants to be determined and subject to 

b 1 + ®2 + * * * +B 7 — 1 (55) 

By letting the components Z 2 » ***, Z 7 satisfy the final conditions 
(53b) , we get 


’ Y 12 

• 

Y 22 ••• Y72 
• • 


' Bx' 
• 


’ m 2 ' 

• 

• 

• 

. yi7 

• • • • • 

• • 

Y 27 * * * Y 77 


• 

. 67 - 


• 

M 7 


Note that z i=3q is related to qi, q 2 » q 3 according to Eq.<3) We solve 
Eqs. (55-56), and eliminate B 7 , to get 

[ Bx B 2 ••• B 6 ]^=[pxo ^ p 10 (N+1) P10 (N+1) r 0 (N+ 1 ) 3 T 
To summarize the above results, we have a step by step description of 
the quasilinearization algorithm: - - 

(a) Obtain nominal function, z^®)(t), and the slewing time t(°)f 
discussed in sections 4.1 and 4.2. 

(b) At the (N+l ) st stage of the iteration (N=l,2, ••*) determine 
whether the components of the control vector are on their 
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or 


bounds (have maximum admissible amplitudes) at any point in 
[0, tf(k)] using the conditions in Eq.(23) .Evaluate the 
matrix A and vector D of Eq.(52) in [0, tf( k )]. 

(c) Solve the linear two-point boundary-value problem, Eq.(52), 
using the method of particular solutions. Obtain the func- 
tion z(N+l)(t). Compute the function u( N+ D(t) from Eqs. 

(49) and (51). 

(d) Compute 

f - « a x{|P i0 ( N+1 >-Pi 0 tN) M r j 0 (,,+1) - r j 0 <N) | . i,j =1,2,3} 

if an< * 1 r jO £1 ; 


? - 


max 


P i0 (N + l).p i0 (N)| |r j0 (H+l).r j0 


<N) 


P j0 (N+1) 


r i t j —1 f 2 , 3 ^ 


! r j0 (N+1) | 
if |P i0 < N+1) |>l and | r j 0 ( N+ D| >i . 

(e) If f<$, a small preselected positive quantity, the optimal 
solution for tf^ k ) is obtained. Then go to (f). Otherwise, 
go to (b) and the iterated solution is used as the nominal 
function for the next stage of the iteration (N^N+l). 

(f) If one of the controls is of a bang-bang type, stop computa- 
tion, and the tf(^) is the minimum time tf*. If not, shorten 
tf (k) to tf( J - f U by a reasonable amount and go to (b). 


3 -18 


5. Numerical Results 


Finally, we apply these methods described in the previous sec- 
tions to the SCOLE slewing motion HI . Fig. 1 shows the configuration 
of the SCOLE. It is composed of a Space Shuttle and a large reflecting 
antenna. The antenna is attached to the Shuttle by a flexible beam. 
Since we only consider the motion of the rigid SCOLE in this paper, 

the flexibility of the beam is ignored. The X, Y, Z axes are the 

y 

spacecraft axes corresponding to roll, pitch and yaw axes, respec- 
tively. The controls considered in this paper include three moments 
about the X, Y, Z axes of the system and two forces applied at the 
center of the reflector in the X, Y directions only. The inertia 
parameters of the SCOLE and the saturation levels of the controls are 
listed in Table 1. 

' Table 1. Inertia Parameters and Limits of Controls 


I H=1, 132 ,508 
I 12=" 7,555 


I 22 ss 7 , 007 , 447 i 33 =7,113,962 

I 23= 115,202 I 3i= 52,293 


SLUG 

-FT 2 


U 1 

u 

xmax 

= 

_ u 

xmxn 

U 2 

u 

ymax 

= 

_ u 

ymm 

U 3 

u 

zmax 

= 

- u 

zmin 

U 4 

^xmax 

= 

“ ^ xmin 

U 5 

^ymax 

= 

“ ^ ymin 


= 10,000 


800 


FT-LB 


LB 


The associated alignment matrix B is 


B 


1 0 0 0 130 

010 -130 0 

001 32.5 18.75 
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We have done some numerical simulations for the following cases: 

(a) A diagonal inertia matrix I is used. The control is assumed 
to be provided only by torques on the Shuttle. No control forces on 
the reflector are assumed. The expected rotation is a rotation about 
one of the three principal axes, through 20 deg., from rest to rest. 
The result is exactly the same as that of the theoretical analysis 
discussed earlier in this paper, i.e., the control torque about the 
slewing axis is of a bang-bang type while the others remain zero. 

(b) Extend the inertia matrix in case (a) to a non-diagonal form. 
The expected rotation is a rotation about one of the three spacecraft 
axes and the rotation angle is 20 deg. 

Fig. 2 and Fig*3 give the control torques for the expected rota- 
tions "X- axis slewing" and "Z-axis slewing”, respectively. Fig. 2 
shows that u x is near of a bang-bang type, while u y , u z are not. The 
non-zero contributions of the u y and u 2 are due to the offset of the 
inertia distribution of the SCOLE configuration (non-diagonal matrix 
I). Similar situations are shown in Fig. 3 where u z is near to a bang- 
bang type and others are not. 

The starting value of t^>f for these slewings (X- and Z- axes) 
are t =12. 5749 sec. and t(°)f=31.5166 sec., respectively, by using 
the method in section 4.2. The minimum time-, t*f, we actually obtained 
are t*f=12.57 sec. and t*f=»31.33 sec., respectively. These results 
indicate that the estimated values for t*f are very accurate. 

Fig. 4 shows the changes of attitude angles O x , O y , O z ( 1-2-3 
Euler angles) for the "X-axis slewing" of Fig. 2. The O x changes from 

y 

zero to 20 deg., but 0 y and 0 Z change very little during the slewing 
and finally approach zero. The non-zero changes in 0 y and O z are also 
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due to the offset of the structural distribution of the SCOLE. 

Fig. 5 shows the attitude changes for the "Z-axis slewing" (Fig. 3) 
Unlike the case in Fig. 4, the Ox changes greatly. This change is due 
to the unsymmetrical moments of inertia about the X-axis and Y-axis. 

Fig. 6 shows the control torques for the "X-axis slewing" with 
a slewing time tf=15.37 sec., which is 2.8 sec. more than the minimum 
time t*f=12 . 57 sec. (Fig. 2). The controls are almost linear functions 
of time (rest-to-rest slewing). u x is less than the saturation level, 
and Uy, u z are near zero. From Fig. 6 and Fig. 2, we see that much more 
control effort (approximate 50%) is saved if we increase the slewing 
time a little. Another feature of using a longer slewing time in the 
computation is that it needs less number (4 times) of iterations for 
convergence than by using a minimum slewing time t*f (12 times). These 
properties suggest that, in practical applications of this problem, 
it is not necessary to seek the minimum time, t*f, and the associated 
extremum controls. It is enough to know approximate values of the t*f 
and the controls. 

(c) Following the case (b), we now add two control force actu- 
actors on the reflector, f x and fy. Figs. 7a, 7b and Fig. 8 show the 
control torques, forces and attitude angles for the "X-axis slewing" 
motion. The slewing time t*f is greatly shortened, t*f=3.988 sec. 
(about one third of the slewing time without the forces f x and fy). 

Figs. 9-14 show the controls and attitude angle changes for the 
"z-axis slewing". For the sake of comparison, we use 3 different tf in 
the computation, tf=27.5 sec., tf=26.1 sec. and t*f=20.0 sec. (minimum 
time; recall that t*f=31.33 sec. without f x and f y ). From Figs. 9-11, 
the control torques approach the bang-bang type, and the maximum 
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amplitude of the control forces increases gradually. From Fig. 12 to 
Fig. 14, we can see the obvious increases in 6 X and &y. This is due 
to the increases in u x , Uy , f x and fy. 

(d) Now we consider a general case. Suppose the SCOLE is in an 
Earth orbit and we need the line of sight to be directed toward the 
center of the Earth. The orbital coordinate system (x, y, z) is shown 
in the following figure. 



ATTITUDE OF THE SCOLE 


Suppose, before the slewing, the Y axis of the spacecraft coin- 
cides with the orbital y axis, and the angle difference between x and 
X (or z and Z ) axes is </=7. 897224212 deg. Thus the initial attitude 
quaternion of the spacecraft is 

qo ( 0 ) = cos (*/2) 

q 2 ( 0 ) = sin (*/2 ) , qi( 0 )=q 3 ( 0)=0 

According to Ref.[l], the unit vector along the line of sight in the 
rigid SCOLE coordinate system is 


A 

r los 


f 0.1112447155 


-0.2410302170 


t 0.9641208678 


The direction cosines of the orbital z axis in the body system at the 
initial time are represented by 

F sinoC 1 


z/b - 


0 


coscC 


Therefore the angle between Rlos an< * Z /B at t ^ e initial time is 

6 ^ 08 ( 0 ) — Rlos * Z /B = 0. 96412086 78 cos 0C — 0 • 1112447155 sincC 


* 20 deg. 

The eigen axis of the expected rotation in the body system is deter- 
mined by 


r LOS x z/B 
r LOS x z/b 


Thus the quaternion for this rotation is 

qo = cos (20^2) 

qi = sin (20^2) , i=l,2,3 

From Eq . ( 5 ) we can get the final attitude quaternion, q(tf). 
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Figs. 15a and 15b show the control torques (without the reflector 
control forces) and attitude angles for this slewing motion. 

The guessed starting value of tf is tf ( 0 )=26 . 3482 sec. and the 
actually converged value of tf is tf*=25.01 sec. The O^os in Fig. 15b 
is the angle between the line of sight and the line of the target 
direction (from the spacecraft to the center of the Earth). 

Conclusion Remarks 

(1) There is a good agreement between the guessed value of tf 
and the value of tf to which the algorithm converges in 
the cases (b) and (d). 

(2) The guessed initial values of the costates here: p(0), r(0) 
are adequate for the algorithm to converge. If the slewing 
time, tf, is sufficiently larger than the minimum time, 
tf*, then, the converged values of p(0) and r(0) are very 
close to the guessed values and less number (4 times) of 
iterations is needed (Fig. 6). The same situation was 
observed in Ref . 8 . 

(3) The control profiles obtained in this paper give us a good 
reference for future use. For example, an extension to the 
minimum time slewing motion of the SCOLE model containing 
both rigid and flexible components is planned. 
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Figure 1. Drawing of the Shuttle/Antenna Configuration. 
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where the aij and b^j (the elements of matrices A and B) are constants 
associated with the inertia parameters of the spacecraft. 

Then the term r T I“lu)IttJ of the Hamiltonian, H, in Eq.(14) has 


tac x. uxui 


r T I -1 tOltO = [ ri r2 r3 ] [ A : B ] iO 
2 2 2 

Rl 1 + R2 ^ 2 +R 3 ^3 +R 4 61) 2 W3+R5 W 3 £*^1 +R 6 ^ 1 £^2 


where R^ are the linear combinations of rj 

R i I 1 

R 2 = A t r ; ] 


R 4 

R 5 = B t r 


The term g(0, r) in Eq.(16) is obtained by 

2Ri R6 r 5 

g((0, r) = - - - r 6 2r 2 r 4 


r 5 r 4 2r 3 
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Appendix II 


Solution of Eg. (44) 


Eq. (44) can be rewritten as 

© = bi © 2 + cf T ^ (44) 

For simplicity, we only consider the solutions for the following boun- 
dary conditions 

9(0H0, e( 0)=0 (II-l) 

0 (t f )=©*, ©(t f )=0 ( II -2 ) 

Suppose aj^O, and let 

b a bi/ai , c = ci/ai 

we can rewrite Eq.(44) as 

e = b e 2 + c T ( ii-3 ) 

where we also suppose c> 0 . 

Since the control for this problem is of a bang-bang type with 
only one switching point, then, for T=+l, 

6 = b 0 2 + c ( II-4 ) 

By integrating Eq.(II-4) and using condition ( II— 1 ) , we can get 

9 = ^(c/b)(e 2bd - 1) ( II-5 ) 

For T =-l. 


© = b 0 2 - c ( II- 6 ) 

By integrating Eq.(II- 6 ) and using condition (II-2), 

0 = ( c/b ) ( 1 - e 2b (3-0* ) ) ( II-7 ) 

By equating Eqs. (II-5) and (II-7), we get ©=© s and ©=©3 at the 
switching point, t=t s . 



log 


2 

1 +e - 2 b&* 


(II-8) 


J" (c/b) (e 2be s - 1 ) 


( II-9 ) 
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By integrating (II-5)/ we get 


V-bc 


cosh“l(e b0 s) , b<0 ; 


( 11 - 10 ) 


[ - sin“l (e“ b0 s)] , b>0 

v'EE 2 


and by integrating Eq.(II-7), we get 

1 , 7C 


f t s + [ - sin”l(eb( 0 *“ 0 s) ) ] , b<0; 

2 


tf - 


1 


( 11 - 11 ) 


t s + cosh‘“l[e^ ) (®*" ,0 s) ] , b>0 


J * be 


For the case 


9 ( 0)*0 

more complicated solutions can be obtained, but are not given here. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

For large ordered systems, typical of mathematical models of proposed 
large space structural systems, a finite amount of time will be required to 
evaluate and process the control signals. These input time delays may cause 
instability in the closed loop systems for which the control laws were pre- 
viously designed without considering the effect of the delay. It is seen 
that such instabilities may result even for delays that are only a small 
fraction of the system's fundamental open- loop period. From this study it 
is also observed that even a small amount of inherent (natural) damping 
in the system can increase the amount of delay that can be tolerated with- 
out the system becoming unstable. The control problem with delayed input 
is also formulated in the discrete time domain and an analysis procedure 
is suggested and could provide the basis for future study. 

Pontryagin ' s maximum principle from optimal control theory has been 
applied to formulate slewing control strategies for a general rigid asym- 
metrical spacecraft system. The slewing motion need not be restricted to 
a single -axis maneuver and the final attitude error after the slew can be 
made arbitrarily small. The control effort and minimum slewing time are 
calculated based on a quasi-linearization algorithm for the resulting 
two-point boundary value problem. Euler's eigenaxis rotation theorem 
is used to obtain a nominal trajectory of the attitude slewing motion 
which is used to determine the initial (guessed) co-state values for 
the algcritnm. Numerical examples based on the rigidized, in-orbit 
model of the SCOLE also include the more general reflector line-of-sight 
slewing maneuvers. 


* 


The guessed initial values for the co- states are sufficient to arrive at 
the (converged) values supplied by the algorithm. The methods used here 
may be implemented for practical control sources which may have more con- 
straints. For example, in the neighborhood of the bang-bang control 
switching points, the jump can be replaced by a linear function of time. 
The control profiles obtained provide a reference for future use. An 
extension to the minimum time slewing motion of the SCOLE model containing 
both rigid and flexible components is planned. 
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